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Abstract 



While radiation-pressure cooling is well known 1-5 , the Brillouin scat- 
tering of light from sound is considered an acousto-optical amplification- 
only process [611. It was suggested that cooling could be possible in 



multi-resonance Brillouin systems [9f|l2| when phonons experience lower 
damping than light [12| . However, this regime was not accessible in tra- 
ditional Brillouin systems [6}]Tq] since backscattering enforces high acous- 
tical frequencies associated with high mechanical damping [t]. Recently, 
forward Brillouin scattering [s] in microcavities 11 has allowed access 



to low-frequency acoustical modes where mechanical dissipation is lower 
than optical dissipation, in accordance with the requirements for cool- 
ing Here we experimentally demonstrate cooling via such a forward 
Brillouin process in a microresonator. We show two regimes of operation 
for the Brillouin process: acoustical amplification as is traditional, but 
also for the first time, a Brillouin cooling regime. Cooling is mediated 
by an optical pump, and scattered hght, that beat and electrostrictively 
attenuate the Brownian motion of the mechanical mode. 

Spontaneous Raman- and Brillouin-scattering are common to almost any me- 
dia. Incident photons are annihilated in these processes, giving rise to scattered 
photons at redder Stokes or bluer anti-Stokes frequencies (see Figure [l]i). These 
scattering events lead to the creation or annihilation of phonons respectively. 
While prior research in optomechanics has used forces including centrifugal ra- 
diation pressure pTIslpIs] and optical gradient force 14 -16 to excite mechanical 



motion, it is only recently that Brillouin scattering and electrostriction ||£-11 



were demonstrated in micro cavities. Further, extensive work in radiation pres- 
sure cooling ^-|5j raises the question of whether it is possible to cool phonon 
modes via a Brillouin process, similar to the scattering process available in 
bulk media [7. For achieving Brillouin anti-Stokes cooling the heating Stokes 
line needs to be eliminated. However, since low acoustic frequencies separate 
the pump line from the Stokes and anti-Stokes lines, filtering out the Stokes 
line against the anti-Stokes line requires a rapid transmission change over an 
extremely small frequency difference, which is not easily available in bulk ma- 
terials. Here we transform the energy flow direction in spontaneous Brillouin 
scattering to enable annihilation of phonons and reduction of the effective acous- 
tic mode temperature. Our technique for breaking the material heating-cooling 
symmetry relies on an ultra-high Q microresonator with an asymmetric spec- 
trum of resonances. Cooling is mediated by a pump optical wave, and a scattered 
optical wave, that beat together and electrostrictively attenuate the Brownian 
motion of the mechanical mode. 



Our experiment (Figure [Tp) is based on a silica microsphere resonator 17 
of Qopticai > 10^- We evanescently couple light into the resonator via a ta- 
pered optical fiber |18|. Light is coupled out from the other side of this taper, 
allowing interrogation of the acoustical mode via the light that is scattered 
from it. This spherical resonator supports the three actors that participate in 
the Brillouin cooling process - two optical resonances (pump and anti-Stokes) 
and one acoustical resonance. The acoustical as well as the optical modes are 
of a whispering-gallery type and circulate in unison with considerable over- 
lap. These three modes are coupled via an interplay of photoelastic scattering 
and optical electrostriction (Figure [l^), which in our case allows cooling of the 
acoustical mode. Experimental interrogation of the acoustical mode is made 
possible through observations of the pump light that it scatters (and Doppler 
shifts) into the anti-Stokes mode. The beat note between the pump light and 
the anti-Stokes light is measured on a photodetector at the output of the taper 
and is a measurement of the acoustical mode. The mechanical deformation of 
the sphere is illustrated with 12 acoustical wavelengths along the circumference 
in Figu re p!p . This illustration corresponds to an experimentally observed (see 
Figures [2j^ acoustical mode of frequency Via = 95 MHz on a 77.5 pm radius 
sphere, which we numerically calculate in the Supplement. In this work, we will 
experimentally cool this mechanical mode. 



While energy considerations give a general explanation for cooling, an ex- 
plicit solution is needed. We therefore also analytically derive the cooling process 
while starting from the acoustical- and optical-wave equations in the Supple- 
ment. The source terms in these equations are taken to be electrostrictive 
pressure for the acoustical wave, and photoelastic index change for the optical 
wave. As expected from energy considerations, the solution shows that light 
applies compressive pressure on the region of the acoustical wave that expands, 
to take away energy from sound. We also confirm that repeating the same cal- 
culation for the Stokes process gives pressure on the acoustical regions while 



2 



they shrink so that energy is imparted to the acoustical wave as expected. The 
classical analysis also shows that higher acoustical and optical quality factors, 
and lower acoustical frequencies, improve the cooling ratio (Eq. S.44). A quan- 
tum analysis 19 provides more details on the feasibility of this system towards 
ground-state cooling 4ppQ| . In what follows we will explain our technique that 



enables anti-Stokes cooling (Figure [T^), without the degradation caused by the 
heating Stokes process. 



The major enablers for Brillouin cooling are (i) the selective resonant en- 
hancement of the anti-Stokes cooling process over the Stokes process, and (ii) 
phase matching for the two optical modes and the acoustic mode interacting 
through the Brillouin process. Selective resonant enhancement of the anti-Stokes 
process (Figure [l^) is made possible due to the non-periodic frequency separa- 
tion [21,22 of the high order optical resonances of the sphere. Here we exploit 
this non-periodic resonance structure for solely selecting the anti-Stokes line. 
The phase matching requirement implies that the anti-Stokes photon should 
carry out the energy (uj) and momentum (M) of the pump photon and the 
thermal phonon. In a spherical resonator, the momentum parameter M relates 
to azimuthal propagation e^^^^~^'^^ around the sphere equator. As expected, 
these conservation considerations are in agreement with our analytically derived 
synchronous solution for the coupled wave equations (?) that describe our system 
(See Supplement). Both phase matching and selective resonant enhancement 
are illustrated in the energy- momentum (uj-M) diagram of Figure [l]3. Since the 
pump optical resonance (cjp, Mp) and the anti-Stokes resonance {uJaS^ Mas) 
are separated by the acoustical resonance parameters (1^^, Ma), as indicated by 
the solid triangles, both phase match and resonant enhancement are satisfied 
for anti-Stokes scattering. On the other hand, the dashed triangle of Figure 
[l]3 illustrates the off-resonantly eliminated Stokes process, where subtraction of 
the same acoustical parameters (1^^, Ma) from the pump resonance (cjp, Mp) 
brings us to a region (a;^, Ms) where no optical resonances exist. 



Experimental measurement of Brillouin cooling is shown in Figure [2| as in- 
dicated by the broadening of the beating signal between pump and anti-Stokes 
lines, as a function of increasing pump power. This beat observed from the 
tapered fiber coupler serves as a measurement of the acoustical mode (See Sup- 
plement). The cooling experiment is performed by positioning the pump laser 
at the lower frequency resonance {ujp in Figure [T]3), and by observing the light 
scattered in the anti-Stokes direction into the uJaS resonance. As in [1^23^ we 
prefer linewidth as a measurement of cooling p^[|5| [2Q||23||25] since the alter- 
native of measuring the integrated power of the beat note might be affected 
by attenuation of the optical signals. Since the power refiectance spectrum is 
proportional to the squared displacement spectrum of the acoustical mode (See 
Supplement), the acoustical linewidth has been extracted from this data and 
is presented in Figure [3^. Here, the linewidth of the 95 MHz acoustical mode 
increases from 7.7 kHz to 118 kHz with increasing pump power, in an open 
laboratory environment at 294 K. We measure the effective mode temperature 
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through its inverse proportionaUty to the hnewidth of the acoustical mode (See 
Supplement). The result indicates an achieved cooling ratio of 15, implying an 
effective mode temperature of 19 K at the maximum observed cooling point. 



The cavity operates in the thermally stable regime 26 . The cooling factor we 
achieved was limited by the optical power available from the pump laser. To 
complete the analysis, we now invert this process by moving the pump to the 
higher frequency resonance. The low frequency optical mode now functions as 
a Stokes resonance for the higher frequency pump. As expected, line narrowing 
is now observed as indicated in Figure [3^, accompanied by growth of the Stokes 
signal (Figure [2j bottom). The Stokes process linewidth and the anti-Stokes 
cooling linewidth converge at low input powers, indicating the mechanical qual- 
ity factor of the mode. This point of convergence provides the 7.7 kHz acoustical 
linewidth calibration at 294 K in Figure |3j so that mode temperature can be 
calculated during the cooling experiment. 

We experimentally confirm that two optical resonances exist that are sepa- 
rated by the mechanical frequency, as required for phase matched cooling. We 
already know that the pump to anti-Stokes beat, measured by a photodetector, 
provides us with a frequency of 95 MHz for the mechanical mode (Figure |4|l, 
right). We then perform the following two independent measurements that 
obtain no information from this electrical measurement. Measuring resonator 
transmission while scanning the pump frequency (Figure [IJd) indeed confirms the 
existence of two optical resonances spaced by 95 MHz to within measurement 
error. We also verify that the relevant resonance structure is indeed asymmet- 
ric relative to the pump by confirming that no resonance appears on the Stokes 
side of pump resonance Op at a spacing of 95 MHz. The optical line frequency 
spacing is also measured by repeating the Stokes experiment 11 but now with 
an optical spectrum analyzer at the output. As expected, the spectrum ana- 
lyzer resolves two optical lines that are separated by 95 MHz (Figure [4^, left) 
in agreement with the measured transmission profile. These three independent 
measurements - beat spectrum, optical transmission, and the optical output 
spectrum - confirm the existence of the resonant triplet that was used in the 
cooling process. 

Here we use forward Brillouin scattering [sjjlT] in a system involving two 
optical modes and one acoustical mode to demonstrate cooling. In comparison, 
previous optical cooling studies have used fluorescence |27 , radiation pressure 
forces [l]-[5], and photothermal pressure [23^ While the end result here com- 



bines a net annihilation of phonons [lj-[5j 20 23 27 in a multi-resonance sys- 
tem [9-12 , no detuning of the pump with respect to the resonance is required 



for achieving cooling. Additionally, the fact that the mechanical mode is of 
a whispering-gallery type here, suggests low dissipation via the support |11 . 
While we selected the 95 MHz mode in this study because of the greater cool- 
ing ratio achievable for such relatively low frequencies (Eq. S.44), mechanical 
modes spanning from tens of MHz 11 to tens of GHz (9j[l0] can be accessed 
with this method. In a broad context, Brillouin scattering belongs to a family of 
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scattering processes that includes Raman- and Rayleigh-scattering (Figure [T]i) . 
Reversing the energy flow direction here, in respect with Brillouin systems such 
[8-11 , raises the question of whether similar inversion is possible with Ra- 
in spite of 



man scattering [28] non-optomechanical cooling of solids 27 29 30 
the higher frequencies involved. 



as 
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Figure 1: Concept: (a) Cooling of an acoustic wave through the interplay between anti-Stokes photoelastic 
scattering and electrostrictive pressure, (b) Selectively enabling Brillouin anti-Stokes cooling (Fig. a), 
while rejecting Stokes heating is achieved by relying on a resonator with asymmetric optical resonance structure 
(right) that matches the energy and momentum of the acoustical mode being cooled (left), (c) Pump photons 
photoelastically scatter from the Brownian acoustical mode, and thermal energy is removed since only the anti- 
Stokes process is resonantly enhanced. Inset: FEM calculation of an acoustical mode with 12 wavelengths around 
the circumference, corresponds to the frequency Qa = 95 MHz on a 77.5 pm radius sphere. The amplitude and 
index perturbation are calculated for 300 K. (d) The family of material- level scattering processes (after 7 ). 
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Figure 2: Experimental observation of cooling and heating of a 95 MHz 
acoustical mode in the microsphere system. The beating spectra between the 
pump and scattered optical signals are presented here, which are a measurement 
of the acoustical mode. Cooling is observed in the broadening of this spec- 
trum as a function of increasing pump power (top). To cool the acoustical mode, 
we pump the low frequency optical resonance (Op) in order to only resonantly 
enhance the anti-Stokes process while suppressing the Stokes process. For the 
heating experiment (bottom), we pump the higher frequency optical resonance 
{Oas) instead, to resonantly enhance the Stokes process. 
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Figure 3: Experimental result: (a) We capture linewidth data from a 
cooling and heating experiment, and calculate the effective acoustical mode tem- 
peratures. As expected, both cooling and heating linewidths converge towards 
the Brownian room temperature linewidth at low optical pump power. Dashed 
line is a guide for the eye. (b) We calibrate the acoustical linewidth as a function 
of temperature by finding the point of convergence of the cooling and heating 
data. The narrowest observed cooling linewidth (top) and the widest observed 
heating linewidth (bottom) are averaged to estimate the 7.7 kHz linewidth for 
the acoustical mode at ambient temperature of 294 K. 



10 



Optical Spectrum Analyzer 
resolves pump and Stokes signals 



Electrical Spectrum Analyzer 
measures beat note 



Stokes 
scattered 
light 



Laser 
pump 



95 MHz spacing 



|j|yi||^m||||ii|lij|^^ mjuM J^n||||||I||||I|uL 



-400 -200 200 400 
Relative optical frequency Aw, MHz 

(a) 



16 



95.1 95.15 
Frequency, MHz 



0.8 



C0.6 



E 0.4 



CO 

^ 0.2 



Op] 



Pump 
resonance 

7 



95 MHz spacing 



A. 




20 40 60 80 

Relative optical frequency Aco, MHz 



100 



(b) 



Figure 4: The participating resonances are analyzed via three indepen- 
dent experimental measurements, confirming that (a) the optical signals 
beating on the photodetector (right) indeed originate from one pump and one 
Stokes optical line as resolved with an optical spectrum analyzer (left), (b) The 
separation between these optical lines corresponds with the transmission spec- 
trum of our device, measured by sweeping the pump wavelength. Absorption 
peaks that overlap the lines resolved in (a, left) are seen. Resonances are named 
according to their role in the cooling experiments. 
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1 Numerical calculation of acoustical mode and 
photoelastic grating 

A question arises on the photoelastic mechanism that turns strain and deforma- 
tion, associated with the acoustical mode, into a traveling index perturbation 
that scatters and Doppler shifts light. Two effects must be considered in this 
calculation - (i) the strain-induced change of refractive index [S37], and (ii) the 
shifting of the silica-air boundary as a function of the deformation [S38]. 

Finite-element simulations (Figure [5| are used to numerically calculate a 
surface mode at the experimentally observed acoustical frequency of 95 MHz. 
We then equate the numerically calculated strain energy of the acoustical mode 
to /c^T, to determine the mode amplitude at room temperature. We have 
determined that our Ma = 12, = 95 MHz acoustical mode has a Brownian 
amplitude Uq = 9.57 fm at room temperature (see section [4|. 

In the second stage of this calculation, the strain and deformation outputs 
of the numerical solver are used to calculate the index change using [S37], [S38]. 
Accordingly, the refractive-index perturbation of the Brownian photoelastic 
grating is estimated to have a maximum amplitude of An/n = 3.4 x 10~^^. 
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2 Analytical derivation of anti-Stokes Brillouin 
cooling 

2.1 Problem setup 

Here we will use coupled- mode theory [S36] to couple two optical resonances and 
one mechanical resonance using the optical wave equation and the acoustical 
wave equation. Our analysis assumes two optical modes with similar polariza- 
tion, which allows for considering only the scalar value of the field. 

Consider two propagating optical fields Ep and Eas within the spherical 
resonator representing 'Pump' and 'anti-Stokes' modes respectively. The fields 
are defined in a spherical coordinate system (r, 6>, (j)) as 

Ep = Ap{t) Ep{r, 0) cos(Mp (jy-uopt) (1) 

EaS = Aas(t) Eas{r, 6) COs{MaS (/) - UJaS t ^ t/j) (2) 

where Ap{t)^ Aas{t) are the real and positive amplitudes of the waves, and 
Ep{r, 6>), Eas{r, 0) are the field distributions in the plane normal to propagation 



Displacement Index An/n 




Uo = 9.57 fm 



max(An/n) 
= 3.4x|0-'o 



Figure 5: Finite-element calculation: The Ma = 12, — 95 MHz mechan- 
ical mode is numerically calculated. The mechanical amplitude is calculated to 
be 9.57 fm as a result of Brownian phonon occupation of the mode at room 
temeprature. The photoelastic grating (An/n) is then numerically calculated 
using the effect of strain. 
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Figure 6: Illustration of spatial mode profiles: The cross- sectional acoustic 
mode profile is represented as 7^(r, 6>), while the profiles for the pump and anti- 
Stokes optical modes can be represented as Ep{r^O) and Eas{r,0). We define 
the cross-sectional area transverse to propagation as S. 



(see Figure [6|. The fields propagate in the azimuthal direction (</)) with integer 
propagation constants of Mp and Mas^ and optical frequencies up and Uas- The 
anti-Stokes field is assigned a relative phase shift ip. Later on, this phase will 
be explicitly calculated in order to satisfy a synchronous solution of the coupled 
wave equations. This phase will turn out to mathematically express the major 
deliverable in this paper as it dictates the energy flow direction from sound to 
light (and vice versa if oOaS is replaced with Stokes frequency). The sinusoidal 
representations are chosen instead of phasors for mathematical convenience. 

The acoustic density wave p is also defined without the limitation of gener- 
ality as 



p = Po^Ap 

= Po^ p(t) n{r, 0) COs{Ma (t^-^at) 



(3) 



where po is the unperturbed density of the material, p{t) is the amplitude of 
the wave, 7^(r, 0) spatial density distribution. Ma is the acoustic propagation 
constant, and ft a is the acoustic frequency. 

It is convenient to choose normalized transverse profiles for the three fields, 
such that 



//. 



■R{r, ef dS = 1 
Ep{r,ef dS = l 



Eas{r,ef dS = l 



(4) 
(5) 
(6) 



where S is the cross-sectional area shown cut away in Figure |6] 
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2.2 Optical wave equation with a photoelastic scattering 
source 

Consider the electromagnetic wave equation 



^2 

Mo 



((n^ + 2n,An)E) (7) 



where we have discarded the higher order terms in = (no + An)^. Here, Uo is 
the refractive index of the material, and An ~ cos(Ma 4^ — t) is the periodic 
modulation photoelastically applied to the index through the acoustic density 
wave. The terms can then be rearranged to obtain 

V^E - ^^Coul-^E = ^oeo-g^ [2noAh Ej (8) 

In our case, E = Ep + Eas is the total electric field. The right hand side of 
this equation is a source term that can be rewritten in the form of photoelastic 
scattering of the optical fields E from the density wave Ap as follows. 



V^E- 



nl d^E iio^ole 



where we have used the relations 



(Ap e) (9) 



Air = 2no An = ( 1^ ) Ap (10) 



.dp, 

7e = Po {de/dp) . (11) 

The electrostrictive constant 7e is a material property governed by the material 
density po and the change of permittivity with change in density (See 7 , chap- 
ter 9). 

The optical wave equation [9] can be approached via coupled-mode theory 
[S36] to determine a synchronous solution for the anti-Stokes wave EaS- For 
this, we only consider the synchronous terms at frequency uJaS- Equation [9] is 
then reduced to 

Equation 12 shows that the coexistence of the Brownian density fluctuation [S] 
and the pump optical mode [l] constitute an optical source term that shines light 
into the optical anti-Stokes mode[2j This mathematical result describes, as ex- 
pected, an experimental reality where the Brownian density fluctuation scatters 
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the pump mode into the anti-Stokes mode while properly Doppler shifting its 
color. 



The L.H.S. of 12 can then be simplified to 



Aas{t) ( - ^ ) Eas{r, e) COs{MaS 4>-tOaSt + ij) 



- ^ 2 uJasKs^t)Eas{r, 0) sm{Mas ^-uJast^^) (13) 

where we make the slowly- varying amplitude approximation A^^g{t) ~ 0. The 
first term on the R.H.S. of equation [13] equals zero as it is the unperturbed 
solution for the anti-Stokes mode. 

We further make the simplification that regions of non-zero optical intensity 
are concentrated at a radial location ^\opt that is almost equal to the sphere 
radius (see Figure [6|. This simplification is acceptable since the transverse 
profile of the optical mode is of the order of 1 pm [S41] while the sphere radius 
is of the order of 100 pm. Equation [13] then becomes 

V'^as - -i ■ -jn^ « --^^A',s{t)Eas{r, 0) sin(M,s 4>-<^ast + ^P) 

(14) 

where the relation c/n = ^opt ^as/^aS has been used. 



Considering the R.H.S. of 12, the product term of Ap and Ep provides the 
simplification 

Ap Ep = p{t)n{r, 0)Ap{t)Ep{r, 0) ■ cos{Ma (l)-Vtat)- cos(Mp (j)-uopt) 



p{t)n{r, 0)Ap{t)Ep{r, 0) • - cos((Mp + 



{Ujp + na)t) (15) 



where in the last step only the term providing for anti-Stokes scattering is 
preserved. In order to satisfy a synchronous solution to [T2j the following two 
relations appear 



Mas = Mp^Ma . 



(16) 
(17) 



This result is expected from energy and momentum conservation considerations. 
Equations [T6] and [Tt] describe a situation where the created anti-Stokes photon 
carries the momentum and energy of the annihilated thermal phonon and pump 
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photon. This result is apphed to the R.H.S. of 12 to give 



(Ap Ep) 



1 ^2 



p{t)n{r, e)Ap{t)Ep{r, 9) ■ -—- cos((Mp + Ma)<j) - (up + Qa)t) 
Po 2 at'^ 



Pof-ole 
2po 



p{t)n{r,e)Ap{t)Ep{r,0) loIs cos{Mas <I> - coas t) . 



(18) 



where the slowly-varying amplitude approximation has been applied. 



Combining the L.H.S. and R.H.S. of 12 (from 



14 



and 



18), multiplying both 



sides by Easix^^)^ integrating over the area we obtain 



^2 Ks(t) MMaS <i)-LOast + ^) / / Eas{r, 0)' dS 



Pof-ole 
2po 



p{t)Ap{t) Lols COs(MaS (/> - W, 



aS t) jj n{r,e)Ep{r,e)Eas{r,e) dS 

(19) 



The normahzation of equation |6] is then apphed, and we reach the result 

A'^sit) Sm{MaS 4> - OJaS t + ip) 

= Ci p{t) cos{Mas 4>-0Jast) [[ 7^(r, 0)Ep{r, 0)Eas{r, 0) dS 



where Ci is the positive pre-factor 

4 Pp Ml, 

Without loss of generality, we assume that the modal overlap integral 



■ Ap{t) . 



^overlap = j j n{r,e)Ep{r,e)Eas{r,e) dS 



(20) 



(21) 



(22) 



is a positive number. A synchronous solution for equation 20 (and the optical 
wave equation 12) is then obtained for the phase angle t/; = 7r/2. The simplified 
form of [20] is then 



^aS'(^) = C'l p{t) ^overlap 



(23) 
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2.3 Acoustic wave equation with electrostrictive pressure 
sources 



The acoustic wave equation can be written as 



Ma 



where the term within parentheses is the velocity of sound, and 

1 



Pst 



-eo7e\Ep ^ Eas\ 



(24) 



(25) 



is the electrostrictive pressure generated by the total electric field (See | 7 , chap- 
ter 9). Similar to the treatment for the optical wave equation, the L.H.S. of the 
wave equation |24] can be simplified to 



Ma 



V^p ^ 2 n„ p'{t) n{r,e) sin(M„ <l> - t) 



(26) 



where some terms have been removed through the use of the slowly- varying am- 
plitude approximation, as well as knowledge of the unperturbed acoustic wave 
equation. 



Considering the R.H.S. of the acoustic wave equation 24 we have 

\Ep + EaS? 

= \Ap{t)Ep{r,e)cOs{Mp <f) - up t) + Aas{t)Eas{r,e)cOs{MaS <j) - UJaS t + 

(27) 

which is simplified to 
\Ep + Eas? ^ Ap{t)Aas{t)Ep{r,e)Eas{r,e) 

(cOs((MaS - Mp) (j) - (WaS " OJp)t + V) 
+ cos((Mas + Mp) - {u)aS + Wp)t + t^)) (28) 
where the high optical frequency terms at 2wp and 2uJaS, as well as the do 



and 17 and preserve the synchronous term from|28]to match equation [26 
1 

eole\Ep + EasV 

at 



terms, have been discarded. Again, we employ the conservation relations [16 

eo-ie\Ep ^ Eas? 



1 

~ ~2 
1 



(29) 



eoieAp(i)Aas(t)Ep{r, 0)Eas{r, 0) cos{{Mas - Mp) - {uoas - ^p)t + V^) 
eo-ieAp{t)Aas{t)Ep{r, 0)Eas{r, 0) cos(M, c^-^at^^) 
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The acoustical-source term (R.H.S. of 24) is then evaluated as the Laplacian 
of the electrostrictive pressure using the spherical coordinate system. 

\/^Pst= (\dr{r^dr)^^^de{smOde)^^^r^dl)pst (30) 

sma' sm 6> / 

Although the common terminology for this term is the 'source', it will be shown 



here to actually attenuate the Brownian density fluctuation. The R.H.S. of 30 



yields several terms including spatial derivatives of the electric field distribution. 



We now reassemble the two sides of the acoustic wave equation ( 24 ) from 26 
and|30j multiply both sides by 7^(r, and integrate over the area S described 
in Figure [gJ Using the normalization for 7^(r, 0) shown in [4j this gives us the 
equation 

2 p'(t) sm{Ma ^-^at) = jj^ {\/%t) ^{r, 0) dS . (31) 

Assuming that the optical modes are confined to the equatorial region as shown 
in Fig ure [6} we have numerically determined that the last term of the Laplacian 
in 30 i.e. (1/r^ sin^ 6>)9^, is dominant in the solution of the R.H.S. integral in 
|3Tj~We therefore obtain 

. ^o^eAAt)AMt)M! .,,,(^^ ff ^P(r,.)K.(r Wr,.) 

2 J Js sin 

(32) 

Next, we apply the approximations r ^ ^opt and 6> 7r/2 for the denominator 
of the integral, since the integral only obtains non-zero values where the optical 
modes are present. We therefore have 

_ C2 Aas(t) cos(M, (^-nt^^) ff 

Vl^^pt Sm{Ma (p - ila t) JJs 

(33) 

where C2 is the positive pre-factor 

C2 = ^^-Ap(t). (34) 

Substituting the overlap integral ^overlap > and the relative phase t/j = 7t/2 
determined from the synchronous solution to the optical wave equation 33, we 
have the simplified form 



// X C2 Agsjt) ^ 

P V^J ^ ^2 -^overlap • \'^^) 

'^opt 
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which indicates that p'{t) is always negative. 

A major outcome of this calculation is now revealed. The synchronous so- 
lution of the two wave equations |9] and [24] implies that when both the pump 
and anti-Stokes optical waves resonate (equations [l][2|, the acoustical resonance 
(equation [3]) is attenuated. 



2.4 Verification of acoustic gain during Stokes scattering 

We will now confirm that the energy flow in [35] is inverted (i.e. goes from light 
to sound) if the Stokes wave is considered instead of the anti-Stokes wave in [2] 
We emphasize that in our cooling experiment the Stokes line was off-resonantly 
attenuated and that this section is only for mathematical confirmation. 

A major difference here is that the Stokes line has frequency ujs lower than 
the pump (as opposed to the anti-Stokes line having higher frequency). Fol- 



lowing the method in section 2.2 for a synchronous solution to the Stokes wave 
defined as 

Es = As{t) Es{r, 6) cos{Ms c^-uost^^) (36) 

we obtain the same result for the phase angle = 7r/2. At this stage there is 
no difference between the solutions for Stokes and anti-Stokes scattering. 

In the acoustic wave equation solved in section |2.3| however, we obtain a 



different result. Searching for a synchronous solution in equation 29 results in 
a change of sign for the ip term in the electrostrictive pressure 



PstU = --eoleAp(t)Aas(t)Ep{r,e)Eas{r,e) cos(M, ^ - t - ^) . (37) 



This is a major difference that eventually results in a sign inversion in the rate- 
of-change of acoustic wave amplitude. 



C2 As(t) 

^opt 



p'{t) ~ H ^overlap • (38) 



As expected, this result indicates that acoustic wave sees gain (or, is heated) 
through Stokes scattering, which is consistent with previous derivations of 
[S44] and [6]. 



2.5 Cooled acoustical mode 

We now incorporate the loss of anti-Stokes photons from the cavity into equation 
23, and we set A'^si^) ~ ^ steady state. 

-^overlap ^aS Aas{t) = (39) 

where 5as is a loss rate defined by the loaded optical quality factor [S35] for the 
anti-Stokes optical mode {5as = '^as/'^Qas)- This provides a solution for the 
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steady state anti- Stokes mode amplitude for a given acoustical density fluctua- 
tion. 

Aas{t) = 9l^^p{t) . (40) 

OaS 



Substituting 40 in 35 reveals that 

fU\ ^1^2 ^overlap /,x 
P W = -Pit) 

/ 7e ^overlap ^ .2] u\ - ^60!!. ... 

= - [ Sc-p^Ml,^^ ■ ^"^^^^ j = - ^ 

where we deflned the R.H.S. with rBriii. being the rate of Brillouin cooling by 
light. 

We now take into account intrinsic mechanical damping and thermal 
fluctuations that enter the equation of motion for the quadrature amplitude p 
as a Langevin noise force term [S43]. The solution of the equation of motion 
reveals that the acoustical density-fluctuation p has a power spectral density of 
Lorentzian shape, with linewidth + rBriii.- As a result, we obtain 

_ (p') _ ^^^^ 



^r.t. {Pit) rBriii. 

where ( ) represents time average, and (p^.t.) Brownian density-fluctuation 

at room temperature T^x.- Here Tgff is the effective temperature of the acous- 
tical mode. Previous studies [S34], 23 , [S42], 1 have also arrived at similar 
expressions relating linewidth to temperature. 

We can substitute the pump fleld-intensity in the cavity (equation 41) 
with input pump power Pinput instead, as 

Pinput = i^Al (43) 

where Qp is the quality-factor of the pump mode, and 77 is the vacuum impedance. 

as 



We then write equation 42 using 41 



^ = JL (44) 

^r.t. 7e ^overlap ^aS D D D P. 

where we used the relationship between acoustical-loss and acoustical quality- 
factor Fa = ^a/Qa- 

This result shows that the cooling due to the anti- Stokes process improves 
with increasing pump optical power Pinput- Increasing optical quality factors 
QaS^ Qpi and the acoustical quality factor Q^, also increases the cooling effi- 
ciency. 
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3 Experimental measurement of the effective acous- 
tical mode temperature 

We experimentally measure the beat note between the pump optical signal and 
the ant i- Stokes scattered signal on a photodetector. We will now discuss how 
this beat allows us to observe the acoustical mode and measure the effective 
temperature Teff. 

We first rewrite [4Q1 as 

/io^oTe ^aS" ^opt ^Dverlap . . 

^aS = . ^ . .2 r Ap-p . (45) 

As a result, the acoustical power spectral density (of p) is also transferred to 
the optical signal Aas around the anti-Stokes frequency oJaS- As the anti-Stokes 
signal now carries the spectral profile of the acoustical mode, its beat note with 
the pump optical signal also has the same spectral shape. 

We then define the reflectance of the pump power Pp ex into the anti- 



Stokes mode Pas oc A^^ as R. Using 45 we get 



R=^^p' (46) 

where the calculation is for powers inside the resonator. The power reflectance R 
is expected to have a Lorentzian spectrum owing to its proportionality with p^ . 
Experimentally, we measure the spectrum of R through the beat note between 
the pump and anti-Stokes signals, and we calculate the linewidth Feff = + 
rBriii. of the cooled density fluctuation. 



We then calculate the effective mode temperature Teff through equation [42 
In our experiment, T^.t. was 294 K, and we measured the room-temperature 
linewidth of the observed 95 MHz mode to be = 7.7 kHz (see Figure 4b). 



4 Brownian amplitude of the acoustical mode 

The thermal equilibrium occupation of phonons (n/c) in an acoustical mode of 
frequency Via is provided by the expression (see [S39] chapter 5) 



(47) 



where /c^ is Boltzmann's constant, and T is the ambient temperature. Ad- 
ditionally, the amplitude of displacement {uq^ as shown in figure |5| from the 
equilibrium position due to the acoustical mode is quantized (see [S39] chapter 
4) based on the expression 



4 n/e 



1 



2 I meff O^a 



(48) 
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where meff is the effective mass of the mode (calculated numerically). Using 
the thermal equilibrium phonon occupancy at room temperature, the above 
equations can be simplified to 



2 



The mode amplitude for the Ma = 12, = 95 MHz mode shown in section [l] 
is calculated to be Uo = 9.57 fm at 300 K. 
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